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ABSTRACT 

The K-band Hubble diagram of Brightest Cluster Galaxies (BCGs) is presented for a large, X-ray 
selected cluster sample extending out to z = 0.8. The controversy over the degree of BCG evolution 
is shown to be due to sample selection, since the BCG luminosity depends upon the cluster environ- 
ment. Selecting only the most X-ray luminous clusters produces a BCG sample which shows, under 
the assumption of an Einstein-de Sitter cosmology, significantly less mass growth than that predicted 
by current semi-analytic galaxy formation models, and significant evidence of any growth only if the 
dominant stellar population of the BCGs formed relatively recently (z < 2.6). 

Subject headings: Galaxies: clusters: general — galaxies: elliptical and lenticular, cD — galaxies: 
evolution — galaxies: formation 



1. INTRODUCTION 

The majority of stars in giant ellipticals found in the 
cores of rich galaxy clusters are old; photometric and spec- 
troscopic studies of cluster galaxies out to z f» 1 suggest a 
formation redshift, Zf, greater than 2, with little variation 
within a cluster, and that secondary bursts of star forma- 
tion account for a small fraction of the stellar mass (e.g. 
Aragon-Salamanca et al. 1993; Ellis et al. 1997; Stanford, 
Eisenhardt & Dickinson 1998; van Dokkum et al. 1998; 
Poggianti et al. 1999). 

However, to understand the process of galaxy forma- 
tion it is necessary to know where the stars were formed 
as well as when. In the traditional view of early-type 
galaxy formation — a "monolithic" collapse at high redshift 
(e.g. Eggen, Lynden-Bell & Sandage 1962; Larson 1969) — 
all the stars were formed in situ, in direct contrast to 
the merger-driven growth of galaxies predicted by semi- 
analytical models for hierarchical cosmologies, such as 
CDM (e.g. Kauffmann & White 1993; Baugh, Cole & 
Frenk 1996). Since the ages of the stars are similar in 
both scenarios, it is the change in mass with look-back 
time that separates the two pictures observationally. The 
current data are inconclusive; for example, the hierarchi- 
cal models are favored by the enhanced merger fraction 
seen in the z = 0.8 cluster MS1054.4-0321 (van Dokkum 
et al. 1999), whilst De Propris et al. (1999) show no ev- 
idence for mass evolution of bright ellipticals in clusters 
out to z » 1. In general it is difficult to follow the evo- 
lutionary history of ellipticals since selection methods can 
seriously bias the samples, c.f. the discussions of progen- 
itor bias in van Dokkum & Franx (1996), the effect of 
preferential selection of the most massive objects at each 
epoch in Kauffmann & Chariot (1998), and the use of color 
selection in Jimenez et al. (1999). 

One approach to minimising such problems is to study 
the evolution of a particular class of ellipticals — brightest 
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cluster galaxies (BCGs) — because of their unique location, 
close to the center of the cluster's potential well. BCGs do 
not appear to be drawn from the same luminosity function 
as other cluster galaxies (e.g. Dressier 1978), which sug- 
gests that they have a distinct formation history. Knowl- 
edge of BCG evolution can therefore provide different con- 
straints on galaxy formation models to studies of the gen- 
eral cluster population. 

The K-band Hubble diagram for BCGs has recently 
been extended to z w 1 by both Collins & Mann (1998; 
hereafter CM98) and Aragon-Salamanca, Baugh & Kauff- 
mann (1998; hereafter ABK98): these observations pro- 
vide the opportunity to measure the luminosity evolu- 
tion of BCGs since evolutionary and pass-band correc- 
tions are insensitive to the recent star-formation history 
of a galaxy at near-IR wavelengths (e.g. Bershady 1995; 
Madau, Pozzetti, & Dickinson 1998). The conclusions 
drawn are contradictory, despite the use of the same cos- 
mology and a common assumption that the stellar pop- 
ulations of BCGs are old and passively evolving; CM98 
assert that the stellar populations of BCGs in the most 
massive clusters have not grown significantly since 
whilst ABK98 argue that their results are in good agree- 
ment with the mass increase of BCGs — by a factor of four 
in an Einstein-de Sitter cosmology — predicted by semi- 
analytical models over the same redshift range. The two 
samples have almost no overlap — CM98 having used an X- 
ray selected cluster catalogue whilst ABK98 used a hetero- 
geneous compilation that was mainly optically selected — 
and it is the aim of the present work to show that the 
results can be reconciled by considering the properties of 
the clusters in the two samples. Section || describes the 
BCG sample used — an extension of that of CM98 — and 
the reduction methods employed, whilst section^ presents 
the results of the analysis and a comparison to those of 
ABK98. Throughout this letter an Einstein-de Sitter cos- 
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mology with Ho = 50 km s 1 Mpc 1 is assumed, and X-ray 
luminosities (L x ) are quoted for the 0.3-3.5 keV pass band. 

2. DATA 
2.1. Sample 

The data presented here comprise K-band observations 
of 76 BCGs. This sample, which incorporates that of 
CM98, spans a redshift range of 0.05 to 0.83, and is drawn 
from the following X-ray selected cluster catalogues: the 
Einstein EMSS (Gioia & Luppino 1994; Nichol et al. 1997; 
Henry 1999), the Southern and Bright SHARC catalogues 
(Burke et al. 1997; Romer et al. 2000), and the ROSAT 
NEP Survey (Henry et al. 1997). X-ray selection is to be 
preferred, since both X-ray luminosity and X-ray temper- 
ature should be more closely related to cluster mass than 
optical richness. 

The X-ray luminosity-redshift coverage of the cluster 
sample is shown by the circles in Figure 1. The additional 
symbols show those clusters from ABK98 with a measured 
X-ray flux or upper limit, except for CI 2155+0334 (also 
known as CI 2157+0347), which has been removed because 
photometric and spectroscopic observations show no evi- 
dence for a cluster (Thimm & Belloni 1994; Oke, Postman 
& Lubin 1998), and CI 0016+16, since it is in both samples. 
The difference in X-ray luminosity coverage at z > 0.5 for 
the two samples is striking; the implications of this are 
discussed in section [| 
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Fig. 1. — Cluster X-ray luminosity as a function of red- 
shift. The circles indicate the BCG sample presented here; 
open for clusters used in CM98 and filled for the new observa- 
tions. Clusters from ABK98 with measured X-ray fluxes are 
shown as plus (+) symbols and the arrow (|) symbol repre- 
sents the 3<r upper limit for CI 1603+4329. The dashed line, at 
L x — 2.3 x 10 44 erg s _1 , shows the luminosity used by CM98 
to separate their sample into high- and low-L^ clusters. 

The K-band BCG observations were made using the IR- 
CAM3 and UFTI cameras on the UKIRT and NSFCAM 
on the IRTF, with some of the data being provided by the 
service programs of both telescopes. IRCAM3 and NSF- 
CAM are 256x256 pixel InSb devices with a field of view 
close to 70" by 70" (the IRCAM3 and NSFCAM pixel 
scales are 0.281 and 0.3 "pixel -1 respectively), and UFTI 
is a 1024x1024 pixel HgCdTe array with a pixel scale of 
0.091 "pixel" 1 , giving a field of view of 92" by 92". The 
observing strategy is the same as presented in CM98: the 
BCGs were imaged using a jitter pattern and separate sky 



exposures were taken for those objects which filled the field 
of view. 

2.2. Reduction 

The data reduction system improves upon that pre- 
sented in CM98, and incorporates elements from the meth- 
ods described in Stanford, Eisenhardt & Dickinson (1995) 
and Hall, Green & Owen (1998). An outline is presented 
below as the method will be fully described in a later pa- 
per. 

The individual frames were masked for bad pixels, dark 
subtracted and divided by the exposure time. A flat 
field image was created by median combination of the 
object images — or separate sky exposures if these were 
available — and applied to the object frames. Masking of 
cosmic ray events was performed on the flattened images 
before they were mosaiced together, which completed the 
processing of those objects with sky exposures. Other- 
wise the mosaic — which is substantially deeper than the 
individual exposures — was used to create an object mask, 
which was then applied to the individual images before 
they were median-combined to form a flat. The flattened 
exposures were then processed as above to create the final 
image. 

Observations of stars from the UKIRT faint standards 
list (Casali & Hawarden 1992) were used to calibrate the 
photometry onto the UKIRT system assuming an extinc- 
tion of 0.088 mag airmass -1 , the median value for K-band 
observations at Mauna Kea. Comparisons of the results 
from repeat observations, both within and between ob- 
serving runs, show that the magnitudes agree to 0.05 mag. 

Aperture magnitudes were measured using a 50 kpc 
diameter aperture and have been corrected for Galac- 
tic absorption using the maps of Schlegel, Finkbeiner, & 
Davis (1998): the correction is small, mostly being less 
than 0.05 mag, but reaching 0.1 mag in several cases. The 
position of the aperture was chosen so as to maximise the 
flux contained within it whilst remaining close to the cen- 
ter of the cluster X-ray emission. Those pixels contami- 
nated by stars and obvious non-cluster galaxies were ex- 
cluded from the calculation, being replaced by values cho- 
sen from regions at the same distance from the aperture 
center. No attempt has been made to remove flux due to 
other cluster galaxies falling within the aperture, and so 
the results are directly comparable to those of ABK98. 

3. RESULTS 

The K-band Hubble diagram for the two BCG samples 
is shown in Figure 2. The lines show model predictions 
calculated using the GISSEL96 code (Bruzual & Char- 
lot 1993), for a solar-metallicity stellar population with 
a Salpeter initial mass function: the solid line indicates 
a no-evolution model for a 10 Gyr old stellar population, 
whereas the other lines are for stellar populations which 
form in an instantaneous burst of star formation at a single 
epoch — Zf = 2 for the dashed line and z/ = 5 for the dot- 
ted line — and then evolve passively. The models have been 
normalised to match the low-redshift, X-ray selected, BCG 
sample of Lynam et al. (1999), following the method used 
in ABK98, assuming a growth curve, d log L/d logr, of 0.7 
for the aperture corrections and a color of R — K = 2.6. 
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Fig. 2. — Magnitude-redshift relation for brightest cluster 
galaxies in the observed K band. Filled and open symbols rep- 
resent those BCGs in high- and low-L^ clusters respectively; 
the division is as in Figure 1. The circles in the top panel 
indicate the sample presented here, whilst the squares are for 
the BCGs in the two z = 1.3 clusters from Rosati et al. (1999), 
where the magnitudes have been measured within 50 kpc diam- 
eter apertures (P. Rosati 1999, private communication). The 
bottom panel shows the sample of ABK98, where the crosses 
are for those clusters without a measured X-ray flux. The no- 
evolution prediction, assuming a 10 Gyr old stellar population, 
is shown by the solid line; passive-evolution models, in which 
the stars form at a single epoch, are shown as dashed (zf = 2) 
and dotted (zf = 5) lines. 



The result remains qualitatively the same as Figure 6 of 
CM98; BCGs in high-L^ clusters form a homogeneous pop- 
ulation which is brighter, and has a smaller scatter, than 
that of \ow-L x clusters. This can be more clearly seen 
in Figure 3, which shows the scatter around the model 
predictions as a function of cluster X-ray luminosity. It 
is this relationship between BCG and cluster properties 
that leads to the contradictory conclusions of CM98 and 
ABK98: out of the eleven z > 0.5 clusters in the lat- 
ter sample, nine have X-ray flux measurements or upper 
limits, with all but two of these having a low X-ray lu- 
minosity (Figure 1). It is unsurprising that these clusters 
are not similar to rich, local clusters, as they were discov- 
ered on the basis of their optical properties (e.g. Castander 
et al. 1994; Holden et al. 1997). The squares in Figure 2 
represent the BCGs in the two z = 1.3 clusters discussed 
by Rosati et al. (1999): the high-L x cluster (solid square) 
was discovered by means of its X-ray emission, whereas the 
low-Z^ cluster (open square) was detected by its galaxy 
population. Although based on only two points, this sug- 



gests that the correlation with environment holds at this 
redshift. 

The semi-analytic models discussed in ABK98 predict a 
factor of ~4-5 increase in the stellar masses of BCGs in 
massive clusters since z = 1, for an Einstcin-dc Sitter uni- 
verse. To test whether the data presented here supports 
this level of evolution, a correlation between redshift and 
the BCG residuals (Am^, e.g. Figure 4) has been sought. 
Passive-evolution models with Zf = 2 and Zf — 5 have 
been used to calculate the residuals — since they provide 
a conservative range for the formation epoch of massive 
cluster ellipticals (e.g. Ellis et al. 1997) — and separate fits 
made to the high- and low-L^ cluster subsamples. Since 
the form of any evolution is unknown a priori, a non- 
parametric rank-order statistic — Kendall's r — was used; it 
also has the advantage that it is insensitive to the choice 
of normalisation adopted for the Bruzual & Chariot mod- 
els. All save one of the fits showed no significant (> 3a) 
evidence for evolution; the exception, at a significance of 
3.6cr, was the high-L^ subsample with zj = 2. To find the 
maximum formation epoch that is still compatible with 
evolution of the high-L^ subsample, Zf was increased from 
2 until the correlation significance dropped below 3cr. Evo- 
lution is found only if the stars formed recently (zf < 2.6). 
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Fig. 3. — Residuals about the model predictions, defined as 
Am k — tobog — ^modei, as a function of cluster X-ray luminos- 
ity. The BCGs presented here are shown as circles, the sample 
of ABK98 is shown as in Figure 1, and the squares represent 
the two clusters from Rosati et al. (1999). The three panels 
are for the models shown in Figure 2: a) no-evolution model 
for a 10 Gyr old stellar population, b) formation at a redshift 
of 2 followed by passive evolution, and c) as for b) but with a 
formation redshift of 5. 
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Fig. 4. — Residuals about the Zf = 5 model for the X-ray 
selected BCG sample. Filled and open symbols indicate BCGs 
in high- and \ow-L x clusters respectively. The lines show the 
expected locus of the residuals for the three models shown in 
Figure 2. 

To quantify the amount of evolution allowed by the 
data, the same parametric form as employed by ABK98 — 
namely M(z) = M(0) x (1 + z) 1 — was used to estimate 
the growth in the stellar mass content of BCGs. Fitting 
for both 7 and M(0) indicates that, in the high-L^ sample 
(which best approximates the cluster selection adopted for 
the semi- analytic models), the typical BCG mass has in- 
creased by a factor of 1.9 ± 0.3 (for Zf = 2) or 1.3 ± 0.2 
(zf = 5) between z = 1 and the present. These growth 
factors are substantially lower than either the factor of 
~4-5 predicted by the semi-analytic models, or the mea- 
sured values of 4.6 (z f = 2) and 3.2 (z f = 5), of ABK98. 
The growth factor can also be estimated by fitting for 7 
alone if one assumes a low-redshift normalisation for the 
model predictions. However, this currently involves apply- 
ing a color-correction to low-redshift optical BCG data, 
which introduces further uncertainty: applying a single 
R — K correction to the X-ray-selected sample of Lynam 
et al. (1999) changes the measured growth factor of the 
high-I^ sample by less than 20%, whilst using the normal- 



isation adopted by ABK98 — based on an optically-selected 
sample — increases the growth factor by 50%. K-band ob- 
servations of the Lynam et al. (1999) sample are being 
obtained to circumvent this problem in future work. 

4. CONCLUSION 

The K-band luminosities of BCGs are correlated with 
their environment: clusters with a high X-ray luminos- 
ity contain BCGs which are brighter, and have a smaller 
scatter, than those BCGs in clusters with a low X-ray lu- 
minosity. The BCG evolution seen by ABK98 has been 
shown to be an artifact of a selection bias in their cluster 
sample; at high redshifts, their clusters are systematically 
less X-ray luminous than their low-redshift sample, and so 
their BCGs are systematically fainter. 

Under the assumption of an Einstein-de Sitter universe, 
non-parametric tests show that the only significant evi- 
dence for BCG mass evolution over the range 0.05 < z < 
0.83 occurs when the dominant stellar population formed 
relatively recently (zj < 2.6). Using the same parametric 
form as ABK98, the masses of BCGs in high-L^ clusters 
are found to have, at most, doubled since 2=1, com- 
pared to the factor of ~ 4 increase predicted, for BCGs in 
massive clusters, by the semi-analytic models discussed by 
ABK98. 
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